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1 Introduction 

In particle collisions, the space-time structure of the hadron emission region can be studied us- 
ing measurements of Bose-Einstein correlations (BEC) between pairs of identical bosons. Since 
the first observation of BEC in proton-antiproton interactions fifty years ago [1 J, a large number 
of measurements have been performed by experiments using different initial states [2,, 3J. At 
the CERN Large Hadron Collider (LHC), BEC were observed for the first time by CMS using 
data at centre-of-mass energies ^ys = 0.9 and 2.36 TeV, collected in 2009 IIH; measurements by 
ALICE at 0.9 TeV were reported in |5J. The present paper reports measurements using data 
taken in 2010 at 0.9 TeV, with a sample increase by a factor 15, and at 7 TeV, for the first time. 
The analysis method is similar to that in m, where more details can be found. In this article the 
results at the two energies are compared and additional studies are performed. 

Constructive interference affects the joint probability for the emission of a pair of identical 
bosons with four-momenta pi and p2- Experimentally, the proximity in phase space between 
final-state particles is quantified by the Lorentz-invariant quantity Q = \/—{pi — pi)^ = 
M-^ — 4m^^ where M is the invariant mass of the two particles, assumed to be pions with 
mass niji. The BEC effect is observed as an enhancement at low Q of the ratio of the Q distribu- 
tions for pairs of identical particles in the same event, to that for pairs of particles in a reference 
sample that by construction is expected to include no BEC effect: 

R(Q) = (dN/dQ)/(dN,ef/dQ). (1) 

The ratio is fitted with the parameterization 

R{Q)=C[l + m{Qr)]-{l + SQ). (2) 

In a static model of particle sources, Q(Qr) is the Fourier transform of the space-time region 
emitting bosons with overlapping wave functions, characterized by an effective size r. The 
parameter A measures the strength of BEC for incoherent boson emission from independent 
sources, S accounts for long-distance correlations, and C is a normalization factor. The correla- 
tion function is often parameterized as an exponential Cl{Qr) = or with a Gaussian form 
^{Q^) = 6^(2')^. Other forms have also been used ([6l and references therein), and several of 
these are mentioned below. In addition a formulation aimed at describing the time evolution 
of the source [|ZJ is considered and compared with the data. 

2 Data and Track Selection 

A detailed description of the CMS detector can be found in fS]. The central feature of the CMS 
apparatus is a superconducting solenoid of 6 m internal diameter, providing an axial magnetic 
field of 3.8 T. The inner tracking system is the most relevant detector for the present analysis. 
It is composed of a pixel detector with three barrel layers at radii between 4.4 and 10.2 cm and 
a silicon strip tracker with 10 barrel detection layers extending outwards to a radius of 1.1 m. 
Each system is completed by two endcaps, extending the acceptance up to a pseudorapidity 
= 2.5. The transverse-momentum (pj) resolution, for IGeV charged particles, is between 
0.7% at // = and 2% at\f]\= 2.5. 

Minimum-bias events were selected by requiring activity in both beam scintillator counters Un- 
charged particles are required to have \tj\ < 2.4 and pj > 200 MeV, ensuring that particles 
emitted from the interaction region cross all three barrel layers of the pixel detector and thus 
have good two-track separation. To achieve a high purity of the primary track selection, the 
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3 Definition of Signal and Reference Samples 



trajectories are required to be recor\structed in fits with more than five degrees of freedom 
(^dof) srid x^/^dof < 5.0. The transverse impact parameter with respect to the collision point 
is required to be less than 0.15 cm. The innermost measured point of the track must be within 
20 cm of the beam axis, in order to reduce contamination from electrons and positrons from 
photon conversions in the detector material and secondary particles from decay of long-lived 
hadrons. 

For this analysis a total of 4.2 million events were selected at ^/s = 0.9 TeV, with 51.5 mil- 
lion tracks passing the selection criteria. At 7 TeV, 2.7 million events with 51.7 million tracks 
were selected from data taken during low-intensity runs. Neither of the two energy samples is 
affected by event pileup. Several minimum-bias Monte Carlo (MC) samples were generated, 
followed by detailed detector simulation based on the GEANT4 package 110) . At 0.9 TeV, the 
MC simulations were generated with several PYTHIA6.4 [llj tunes (D6T, DW, PerugiaO, Zl 
and Z2 [fMH ). At 7 TeV, the simulations use PYTHIA6.4 tunes (ProPtO, PerugiaO, Zl and Z2) 
and PYTHIAS. 1 lllSl . 



All pairs of same-sign charged particles with Q between 0.02 and 2 GeV are used for the mea- 
surement. The lower limit is chosen to avoid cases of tracks that are duplicated or not well 
separated, while the upper limit extends far enough beyond the signal region (confined to 
Q < 0.4 GeV) to allow verification of a good match between signal and reference samples. 
The Q resolution in the signal region is better than 10 MeV. Coulomb interactions between 
charged particles modify their relative momentum distribution. This effect, which differs for 
pairs with same charge (repulsion) and opposite charge (attraction), is corrected using Gamow 
factors [16J. As discussed in [4J, the reference sample in the denominator of Eq. ^ can be de- 
fined in several ways: opposite-charge pairs; opposite-hemisphere pairs, where particles are paired 
after inverting the three-momentum of one of them, this procedure being applied to pairs with 
same and opposite charges; rotated particles, where pairs are constructed by inverting the x 
and y components of the three-momentum of one of the two particles; pairs from mixed events. 
In the case of pairs from mixed events, particles from different events are combined with the 
following methods: i) events are mixed at random; ii) events with similar charged-particle mul- 
tiplicities in the same rj regions are selected; iii) events with an invariant mass of all charged 
particles similar to that of the signal are used to form the pairs. In this paper, we use the sam- 
ple obtained by pairing same-sign charged particles from different events that have similar 
charged-particle multiplicities in the same t] regions. This method avoids the possible effect of 
remaining correlations [17J between particles within the same event. The r.m.s. spread of the 
results obtained from the different samples is taken as a conservative systematic uncertainty. In 
dJ an additional, "combined" reference sample was obtained by summing the Q distributions 
of the seven corresponding reference samples. It has been checked that the results obtained 
with the "combined" reference sample are compatible within errors with those presented here. 

In order to reduce possible biases in the construction of the reference sample, a double ratio TZ 
is defined. 



where the subscripts "MC" and "MC, ref " refer to the corresponding distributions from the 
simulated events, generated without BEC effects. 
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4 Determination of Bose-Einstein Correlation Parameters 

Figure [T] shows the distributions of the double ratio TZ for Q > 0.02 GeV and both centre-of- 
mass energies, computed using the tune Z2 of the PYTHIA6.422 simulation, which best de- 
scribes the measured track distributions (in particular the charged-particle multiplicity). The 
shapes fitted with the exponential parameterization Q(Qr) = e^^^ in Eq. (j2| are superimposed 
and the results of the fits are given in Table [l] The values of the two parameters r and A are basi- 
cally related to different features of the distributions: the width of the peak at small Q to r and 
the height to A; they are, however, strongly correlated, with correlation coefficients of about 
86%. The fit quality is poor, as can be seen from the values of ;\;^/Ndof- Gaussian parameteri- 
zations, which are used by some experiments, provide values of X^/^dof larger than 9, which 
confirms the observation in IH that an exponential parameterization is preferred. Compared 
to an exponential shape, alternative functions, as defined in 1181 1191 , and the Levy parameteri- 
zation, Q(Qr) = 6^(6'')" |20|, yield fits of only slightly better quality. However, the latter better 
describes the data at low Q, with values of a < 1 as in |]4|. As a cross-check of the stability of the 
measurement of the width of the peak at small Q and of the fact that it does not depend on the 
fit quality, the average values (first moment) of the Q(Qr) distributions over the same interval 
in Q are found to be consistent for the different functions. More discussion on the shape of 
the TZ distribution and on the fit quality can be found at the end of this section. As discussed 




Q (GeV) Q (GeV) 



Figure 1: Distribution, for Q > 0.02 GeV, of the double ratio TZ defined in Eq. (|3jl, for data at 
a/s = 0.9 (left) and 7 TeV (right). The reference sample is obtained from same-sign charged 
particles from mixed events with similar multiplicities, and the MC simulation is PYTHIA6.4 
tune 72. The lines are the results of fits using the exponential parametrization for Q(Qr), with 
the values of the parameters given in the text. The error bars on the data points are statistical 
only. 



Table 1: Results of of fits using the exponential parametrization for Q(Qr) to the double ratio 
7?., at i/s = 0.9 and 7 TeV. The reference sample is obtained from same-sign charged particles 
from mixed events with similar multiplicities, and the MC simulation is PYTHIA6.4 tune 72. 
Errors are statistical only. 





AVNdof 


C 


A 


r (fm) 


5 (10-^ GeV-^) 


0.9 TeV 


2.5 


0.965 ± 0.001 


0.616 ±0.011 


1.56 ± 0.02 


2.8 ±0.1 


7 TeV 


3.8 


0.971 ± 0.001 


0.618 ± 0.009 


1.89 ± 0.02 


2.2 ±0.1 



in Q, the main experimental uncertainty is due to the choice of the reference sample. In ad- 
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dition, we consider here the systematic uncertainty due to the choice of the MC sample; it is 
obtained from the r.m.s. spread of results obtained using the various MC simulations listed in 
Sect.|2] The results at 7 TeV show a larger dependence on the MC choice than at 0.9 TeV. Finally, 
the uncertainty related to the Coulomb corrections is taken to be ±15%, as determined in [4J, 
which covers the spread from the different parameterizations 121]. The Coulomb corrections 
affect the signal mainly at very low Q, leading to an uncertainty of ±2.8% on A and ±0.8% on 
r. Contributions to the systematic uncertainties are reported in Table |2| and the total systematic 
uncertainties are obtained from their quadratic sum. It was checked that reducing the fit range 
to 0.04 < Q < 2 GeV, thus excluding the first two points at low Q in Fig. [ij gives consistent 
results within errors. The BEC parameters are thus measured to be 

Table 2: Systematic uncertainties on the parameters A and r at two different centre-of-mass 
energies. 





0.9 TeV 


7 TeV 




A 


r (fm) 


A 


r (fm) 


Choice of the reference sample 


0.017 


0.11 


0.015 


0.10 


Choice of MC dataset 


0.009 


0.05 


0.032 


0.16 


Effect of Coulomb corrections 


0.017 


0.01 


0.017 


0.02 


Total 


0.026 


0.12 


0.039 


0.19 



r = 1.56 ± 0.02 (stat.) ± 0.12 (syst.)fm 
A = 0.616 ± 0.011 (stat.) ± 0.026 (syst.) 

at a/s = 0.9 TeV and 

r = 1.89 ±0.02 (stat.) ± 0.19 (syst.) fm 
A = 0.618 ± 0.009 (stat.) ± 0.039 (syst.) 

aty/s = 7 TeV. 

As will be shown below, the increase of r is related to the different average charged-particle 
multiplicities at the two energies, while the value of the A parameter is stable. The BEC sig- 
nal is studied as a function of the charged-particle multiplicity in the event, N^h, as in [4J, and 
of the pair average transverse momentum kj, defined as half of the absolute vector sum of 
the two transverse momenta, kj = |Pyi ± Pt,2|/2. A dependence on kj has been observed at 
the Tevatron [22J and at RHIC II23II , where it is associated with the system collective expan- 
sion. Figure |2] shows the double ratio 7?. as a function of Q for different values of N^h and 
kj. The kj dependence of the r and A parameters for three intervals of multiplicity, obtained 
with the exponential parameterization, is shown in Fig. |3]and given in Table|]forthe0.9 and 
7 TeV data. The effective radius r is observed to increase with multiplicity, for all reference 
samples and for all MC models and tunes, in agreement with previous results. At low mul- 
tiplicity, r is approximately independent of kj, and it decreases with kj as N^h increases. The 
A parameter decreases with increasing multiplicity and kj. The systematic uncertainties are 
estimated to be the same as for the overall measurements (4.2% and 6.3% for A, 7.9% and 
10.1% for r, at 0.9 and 7 TeV, respectively). It should be noted that these uncertainties are 
point-to-point correlated, since the effects of the choice of the various reference samples and 
of the various MC simulations are very similar for the different subsamples. The 0.9 TeV re- 
sults agree with those of ALICE [5]. Figure |4] presents the distribution of the parameter r 
as a function of N^h for both centre-of-mass energies. The measurements are consistent, in- 
dicating that the difference between the values of r obtained for the two global samples are 
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Figure 2: Distributions of the double ratio as a function of Q, for three intervals in kj and 
three intervals in charged-particle multiplicity in the event, N^^, for ^/s = 7 TeV. The lines 
are the results of fits using the exponential parametrization for Q(Qr), with the values of the 
parameters given in Table |3] The error bars are statistical only. 



accounted for by the different average charged-particle multiplicities, which are 12.1 and 19.2 
for the 0.9 and 7 TeV cases, respectively. This trend is consistent with the result of a simi- 
lar comparison between data at 0.9 and 2.36 TeV [4J. The multiplicity dependence is fitted 
as r(Nch) = a ■ N^^^ [24J, giving a = 0.597 ± 0.009 (stat.) ±0.047 (syst.) fm at 0.9 TeV and 
a = 0.612 ± 0.007 (stat.) ± 0.063 (syst.) fm at 7 TeV 

As was noted above and can be deduced from the X^/^do( iri Table [l| none of the quoted 
functions is able to provide a good description of the TZ distributions. This is due to an anti- 
correlation effect between same-sign charged particles for Q values just above the signal region 
(dip with TZ < 1), as shown in Fig.|5j This anticorrelation is observed in the double ratio at both 
energies with any choice of reference sample and MC simulation. It shows little sensitivity to 



6 



4 Determination of Bose-Einstein Correlation Parameters 



E 
























▼ 






CMS - ^ 


\/s = 0.9 TeV 


E 

2 5 








CMS - ^ 


4s = 7 TeV 


2 






T 






2 


■ 




T 






1.5 
1 


■ 

• 




■ 
• 




T 
■ 

• 


1 .5 
1 


• 




■ 

• 




▼ 

■ 

• 


0.5 



~ , 1 , 


. . . 1 , , 


1 


, , 1 




0.5 





, , 1 , , 


1 






0.2 


0.3 


0.4 


0.5 


0.6 0.7 


0.2 


0.3 


0.4 


0.5 


0.6 0.7 

IV-[- V J 


























1 


- • 






CMS - ^ 


>is = 0.9 TeV 


1 


L • 






CMS-^ 


= 7 TeV 




T 












: i 










0.8 


T 
■ 




• 






0.8 






• 






0.6 






■ 

T 




• 

■ 


0.6 






f 




• 

■ 


0.4 


- • 


^ch range: (2 - 


9) 




T 


0.4 


- • 


^oh range: (2 ■ 


9) 




T 


0.2 


■ 


N^i, range: (10 


-24) 






0.2 


■ 


N^^ range: (10 


-24) 









▼ 


N^^ range: (25 


-80) 









T 


N^i, range: (25 


-80) 






0.2 


0.3 


0.4 


0.5 


0.6 0.7 

(GeV) 


0.2 


0.3 


0.4 


0.5 


0.6 0.7 

(GeV) 



Figure 3: Values of the parameters r (top) and A (bottom), as a function of kj in three intervals of 
charged-particle multiplicity in the event, Nch/ for ^Js = 0.9 (left) and 7 TeV (right). The points 
are presented at the position corresponding to the mean value of kj in the considered interval 
of Nch- The error bars are statistical only (in some cases they are smaller than the marker size). 
The systematic uncertainties are discussed in the text. 

Table 3: Results of fits using the exponential parametrization for n(Qr) to the double ratios TZ, 
for three intervals in kj and three intervals in charged-particle multiplicity in the event, N^-h/ 
for v^s = 0.9 and 7 TeV. The errors are statistical only. The systematic uncertainties, which 
are point-to-point correlated, are estimated from the relative uncertainties affecting the overall 
measurements (see text). 
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kj, while it decreases with increasing charged-particle multiplicity in the event, as shown in 
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Figure 4: Parameter r as a function of the charged-particle multiplicity in the event, N^h, for 
^ys = 0.9 and 7 TeV. The dotted and solid lines represent the results of the fits described in 
the text to the 0.9 and 7 TeV data, respectively. The inner error bars represent the statistical 
uncertainties and the outer error bars the statistical and systematic uncertainties, added in 
quadrature. The systematic uncertainties are dominating and are point-to-point correlated. 
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Figure 5: Detail of the distribution of the double ratio TZ for ^/s = 0.9 (left) and 7 TeV (right). 
The dotted lines correspond to fits with Eq. Q, and the solid lines to exponential fits. Note the 
enlarged scale on the y axis. The error bars are statistical only. 



Fig. |6]for the 7 TeV data. This detailed observation is made possible by the large data samples 
studied here, and constitutes the first evidence of this effect at the LHC. Such a structure was 
observed in e^e^ collisions at LEP II25I . The presence of a region of anticorrelation between 
same-sign charged particles has been explained in [7]. In this reference, a parameterization for 
R(Q) has been proposed, aimed at describing the time evolution of the source: 

i?(Q) = C [l + A(cos[(roQ)2 + tan(a7r/4)(Qr)'^]e-(e'-)")] -{l + SQ). (4) 

Here vq is related to the proper time of the onset of particle production, and r enters in both 
the exponential and the oscillation factors. As visible in Fig. |5j fits with Eq. Q are of good 
quality, much better than with an exponential function, and the X'^/^doi is 1-1 both for the 0.9 
and 7 TeV data samples. The depth of the dip in the anticorrelation region is measured as the 
difference A between the baseline curve defined as C • (1 + ^^Q) and the value of R defined 
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Figure 6: Detail of the distribution of the double ratio TZ for = 7 TeV using different intervals 
of charged-particle multiplicity in the event (Nch)- The lines are fits to the data with Eq. Q. The 
error bars are statistical only. 



by Eq. at its minimum. Results are shown in Fig. [7| The depths are found to decrease 
with consistently for the two centre-of-mass energies. The systematic errors have been 
computed from the r.m.s. spread of the results obtained with the various reference samples 
and MC simulations. It has been checked that these results are robust: when the fitting range 
is extended to Q = 5 GeV the results are consistent within errors and the trend is similar. 
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Figure 7: Depth A of the dip in the anticorrelation region, as a function of the charged-particle 
multiplicity in the event, for \/s = 0.9 and 7 TeV. The inner error bars represent the statistical 
errors and the outer error bars the statistical and systematic errors, added in quadrature. The 
systematic uncertainties are dominating and are point-to-point correlated. 
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5 Conclusions 

Bose-Einstein correlations have been measured using data collected with the CMS experiment 
in proton-proton collisions at the LHC, with centre-of-mass energies of 0.9 and 7 TeV. The sig- 
nal is observed as an enhancement of pairs of same-sign charged particles with small relative 
momentum. The parameters are obtained from fits using the exponential parametrization for 
Q(Qr) to the distribution of Q. In agreement with previous results, an increase of the effec- 
tive emission radius r with charged-particle multiplicity in the event is observed, which ac- 
counts for the increase of r from 0.9 to 7 TeV. The parameter r is nearly independent of the 
average transverse momentum of the pair of particles, and decreases with ki in events with 
large charged-particle multiplicity. Anticorrelations between same-sign charged particles are 
observed for Q values above the signal region as previously reported with LEP data. The anti- 
correlation effects decrease with increasing charged-particle multiplicity in the event. 
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